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Abstract 

The solutions for the imaginary susceptibility of the Raman field transition 
with arbitrary relaxation rates and field strengths are examined for three different 
sets of relaxation rates. These rates correspond to (1) Far infrared (FIR) Raman 
lasers in the diabatic collision regime without consideration of coupled population 
decay in a closed system, (2) Raman FIR lasers in the diabatic collision regime with 
coupled population conserving decay, and (3) IR Raman gain in cesium vapor. The 
model is further expanded to include Doppler broadening and used to predict the 
peak gain as a function of detuning for a frequency doubled alexandrite 
laser-pumped cesium vapor gain cell. 
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Introduction 

In this paper we consider the gain in a coherently pumped three level 
system by examining solutions of the density matrix calculation with arbitrary 
relaxation rates, field strengths, and detuningsi 

The problem of a three level laser was first examined quantum mechanically 

by Javan.t 2 ! However, his solution is valid only for a weak stimulated field, or for 

equal pump and stimulated field detuning in the strong field case. Many other 

authors^’^’**’^ have examined the three level system with various approximations, 
most notably Panock and Temkin^l whose solution is applicable for arbitrarily strong 
pump and stimulated fields, but assumes all relaxation rates to be equal. These 
solutions are necessary for the description of resonance enhanced Raman techniques 
for generating tunable radiation. Tunable radiation has been obtained using tunable 
pump lasers and various gases in the far infrared and near infrared regions of the 
spectrum. The far infrared systems utilized C0 2 TEA lasers as tunable pumps and 
molecular gases such as CH S F. The infrared systems were based around tunable dye 
laser pumps and atomic vapors such as the alkali metals. The recent development of 
the tunable alexandrite laser and efficient frequency doubling and tripling crystals 
has made possible solid state laser pumped alkali metal vapor infrared systems. 

In Section I we present the model used to approximate the gain in a 

coherently pumped three level system and the basis of the calculation of the 
susceptibilities. In Section II we examine the effects of different relaxation rates on 
the homogeneous line shape including how the various contributing terms in the 
imaginary stimulated susceptibility saturate with stimulated field intensity. 

In Section III the model is expanded to include Doppler broadening. We then 
briefly examine some of the properties of the Doppler broadened line shape and show 
that although the results at low pump power and detuning within the Doppler width 
agree with the previous calculation,^’**! the line shape for large pump fields and 
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detunings outside of the Doppler width does not reduce to the homogeneous case. We 
then apply the inhomogeneous model to a doubled alexandrite pumped cesium vapor 
laser in order to study stimulated infrared Raman emission. Much of the present 
theory of stimulated Raman emission in non-linear optics is based on calculations of 
from second order perturbation theory.^’ * ^ ^ This approach neglects many 
important effects of the complicated population dynamics associated with the 
relaxation processes of any multi-level quantum system. It is also shown that the 
Doppler broadened prediction of stimulated gain with tuning differs significantly 
from the homogeneous prediction. 
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I. Theory 

We have modeled coherently pumped atomic and molecular lasers of a 
three-level system in an "inverted V" configuration, shown in Fig. 1. The complex 
susceptibilities were obtained^ using the density matrix formalism while neglecting 
all non-resonant terms in the interaction picture. This treatment differs from 

others^’^’^] in that both arbitrary longitudinal and transverse rates as well as 
arbitrary field strengths are allowed. The system is assumed to be closed so that 
there is no net loss of probability. 

The zero field rate of decay of the average probability in any state will 

depend upon the configuration of the entire system and not only upon the value of 
the probability in that particular state. For instance, if the system in Figure 1 was 
prepared in a configuration in which the ensemble averaged probability of finding 
the system in state |2> was equal to its thermal equilibrium value while that of state 
|3> were greater than its thermal equilibrium value, state |2> will not remain in 

equilibrium, since one channel of decay for state |3> is through state |2>. Therefore, 

there will be an increase in the probability of finding the system in state |2>. It is 
clear then that any one state willl retain its equilibrium probability only if the 
entire system is in equilibrium. Rates of decrease or increase of probability in the 
zero field case can only be defined between level pairs so that when fields are 
present the time dependence of the diagonal elements are governed by 

-i - 1 1 

Pll * ft (^31^*13 ' ^13^31^ ^^) " P 33 ) T ^21^11 P22 ) (1) 

'* A 1 *3 

-i ■* 1 1 
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fts = ft |^%lPl3 ' ^13^31^^23^32 "^32^23^ ^ M * ^“^33* ^32^22^" "-p" ^33 ' ^31^11^ 
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where the fy’s are equilibrium occupation ratios (ratio of i th to j th equilibrium 
population), and T r T 2 , and T 3 are the spectroscopically measured relaxation rates of 
the specified level pairs. Notice that trp = 0 is satisfied identically, showing 
probability conservation. 

Rewriting Eqs. (1-3) in terms of the population differences A 13 = p n - p 33 
and A 32 = p 33 -p 2J and Trp = 1, we have 

-1 2 * 

A 1S = £“ Im [^23°31 Es 1 + Im ^31 0 13 E p 1 

- k l [ ( A 13 - O + k 2 [(^32- A 32 )] (4) 

"1 2 

A 32 “ ~~f~ f^31 ” °13 E p 1 + ^ m tl i 23 <7 32 E * 1 + ^3^ A 13 * A 13^ " ^4^ A 32' A 32^ (-*) 


where are the density matrix elements in the interaction picture. E s and E p are 
the complex amplitudes of the signal and pump fields and the constants k r k 2 , k 3 , 
and A 32 are defined in terms of T 2 , T 2 , T 3 and f y by the following relations. 
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The off diagonal elements, on the other hand, give the average amplitude mixing 
between level pairs. The zero field decay of these elements can be approximated to 
first order by a simple exponential decay without contradiction J 
The resulting susceptibilities are given by: 
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where 8 g and 8 are the signal and pump detunings; the T’s are shown in Fig. 1, 

^32^« ^13 ‘ E p 

the a’s, y’s and A’s are defined in Ref. [1], and 8 g = — — » p = — • 

To explore the effects of different relaxation rates, we compare the 
susceptibilities of three different sets of rates at various detunings, pump intensities, 
and signal intensities. The first set of rates correspond to previous results first 
obtained by Panock and Temkin^l in which the decay of each state is uncoupled 
from the other states and all rates are assumed equal. The next set also assumes 
equal rates. However, since the rates are the transition rates between level pairs, 
the pump and signal inversion decays become coupled and the effective decay rates 
are no longer equal. This set of rates is meant to approximate a homogeneously 
broadened molecular gas laser in the infrared or far infrared where energy spacings 
are small compared with the collision energy. In this diabatic collision regime the 
level pair and polarization decays are collision dominated, and thus, all decay rates 
are equal to the collision rate. 

The third system to be modeled is a cesium vapor Raman laser which is 
pumped from the 6S^ 2 level to the 8P^ 2 level and emission is to the 8S x / 2 level. 
The vapor is at low pressure so collisional effects are negligible and decay is 
dominated by spontaneous emission. 

Since the Doppler width of cesium at temperatures above the vaporization 
temperature is much greater than the homogeneous line width, the susceptibilities 
must be Doppler integrated to obtain the actual gain of the system. However, it is 
instructive to first examine the effects of the different relaxation rates on the 
homogeneous line shape, before including inhomogeneous effects to more accurately 
predict the tuning curve of a cesium vapor Raman cell as a tunable infrared source. 
The presence of the 6P X/ , 2 , 7S^ 2 and 7P^ 2 levels complicate matters slightly since the 
system actually has six levels. However, neglecting the spontaneous emission fields 
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and assuming the only fields present are the fields nearly resonant with the pump 
and laser transitions and that all non-resonant contributions of these fields are 
negligible, then there are no field dependent perturbations connected to the 
TSj/j, or TP 1 j 2 levels. The energy levels of the cesium system are shown in Fig. 2. 
The equations of motion of the diagonal elements are as follows: 
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• ^ ^ 
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Where the labels 4, 5, and 6 denote the 7P^ 2 , 7S x y 2 , and 6P 1 y 2 labels, respectively, and 
the zero field equilibrium occupation of the upper levels are negligible. In steady 
state, Eqs. (16-21) can be rewritten, using trp = 1, as two equations in terms of A 13 
and A 32 exactly as in the three-level case except the effective rates k x , k 2 , k 3 , and k 4 
will be different. The dephasing rates in this radiative limit will be taken to be the 
sum of the spontaneous emission rates of each decay channel for both levels 
constituting the level pair in question. The dcphasing rates are given by 
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The spontaneous emission rates are calculated using 


1 
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€hc 3 


(25) 


with the dipole moments obtained from Ref. [13]. The spontaneous emission rates, as 
well as the effective transverse and longitudinal relaxation rates are given in 
Table 1. 


II. Homogeneous Broadening Results 

In this section we compare the imaginary part of the susceptibility at the 
Raman field frequency, which is responsible for gain or attenuation, for different 

sets of rates. It will be shown that many of the qualitative features of the previous 

results, corresponding to the first set of rates, will be changed by the more involved 
relaxation process occurring in the second two systems. The small signal or 

unsaturated case was first examined for equal uncoupled decay rates by Javan^l and 
later Panock and Temkin^ as a special case of their work. It was shown that for 
off resonance pumping (6 p >Av), the susceptibility exhibited two peaks of equal height, 
one near 8 g = 0 and the other near 6 g = 8 p , but displaced by the Rabi shift of the 

energy levels. When the laser field is turned on the symmetry is broken and the 8 g 

24 0 peak saturates faster than the 8=8 peak until there is only the peak which 

B P 

shifts towards the center (Fig. 3(a-d)). Physically, this is because the quantum 
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efficiency of a double quantum transition is one. Or in semiclassical language, it is 
because the single photon, 8 g — 0, contribution is proportional to p u - p 22 , while the 
two photon, 8 g * 8 p is proportional to p u - p 22 . The effect of the saturating field 
to first order is to induce transition from p 33 to p 22 , which reduces the single photon 
inversion, but for any reasonable pump intensities the ground state will have a much 
greater population than either excited state so that the two photon inversion saturates 
much slower than the single photon. 

When the susceptibility was decomposed into single and double photon 
contributions, it was seen that the one photon part was absorptive at the gain peak 
of the two photon part, and similarly, the two photon part was absorptive at the one 
photon peak. The single photon contribution had a much higher peak than the two 
photon contribution, but the two photon absorption at the one photon peak was much 
greater than the one photon absorption at the two photon peak, such that when 
plotted together the 8 g - 0 and the 8 f * 8 p peaks were of equal height (Fig. 4). It 
can also be seen from the decomposed susceptibility that the decrease in gain with 
increasing laser field was due not only to saturation but also to an increasingly more 
effective competition between gain and absorption. It can be seen in Fig. 5 that 
initially only the single photon peak saturates until it is approximately equal to the 
two photon absorption. At this point, the one photon absorption begins to grow 
while the two photon gain then begins to saturate (Fig. 6). 

In addition, the one and two photon contributions nearly cancel, leaving a 
single peak in the total susceptibility between the original 6 g s 0 and 8 g 3 8 p peaks, 
approximately an order of magnitude less than the peak of the single photon 
contribution. The absorptive part of the two photon contribution remains 
approximately the same value as its 0 g = 0 value and only now begins to saturate as 
both contributions saturate uniformly and the resultant peak in the susceptibility 
shifts towards 8 g = 0 (Fig. 7). 
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The on resonance case is harder to interpret in terms of a photon picture, 
especially at powers at which the Rabi splitting becomes significant. At low pump 
powers the one and two photon contributions both have peaks at 8 g = 0, with the 
two photon contribution having absorptive side bands (Fig. 8). At higher powers, the 
single photon term surprisingly becomes a split absorption peak at 8 g = 0, with gain 
on the side bands, while the double photon term has a split gain peak with 
absorption on the side bands (Fig. 9a). The sum of these parts give the familiar 
Rabi split line shape (Fig. 9b). With increasing signal field both parts saturate and 
the splitting and sidebands disappear and the one and two photon contributions are 
respectively single absorption and gain peaks at 8 g = 0. This causes the total 

susceptibility to decrease in height and the splitting is lost (Fig. 10). The scenario 
described above applies only to the equal uncoupled rate case. 

When the rates are unequal, the resonance behavior is qualitatively similar 
but there are significant differences in the off resonance behavior. This is to be 
expected, since the denominators of the population differences A 13 and A 32 contain the 

population relaxation rates k lf k 2 , k 3 , and k 4 , as well as T lt t 2 , and r 3 in the 

different parts, which alter what might be called the effective saturation intensities 
of the different parts. Even for 8 g = 0 there will be differences because of a 

bottlenecking effect which occurs when one of the population relaxation rates is less 
than the others and thus cannot relax fast enough to keep the occupation of that 

level down to the valve it would have in the equal rate case. As a result the 

balance of absorption and gain required to make the 6 g - 0 and the 8 g = 8 p peaks 
equal at 6 g = 0 will not be present. 

In the FIR system with coupled rates, the gain of the one photon term is 
less than the absorption of the two photon term, and therefore, the resulting 

susceptibility has net absorption not gain at 8 g = 0 (Fig. 11). Then, as the signal 
field is turned on, the peak of the one photon term increases to cancel the two 


- 12 - 


photon absorption instead of saturating, as would be expected (Fig. 12). 

In the cesium system the peak of the two photon term is much greater than 
the one photon peak or the two photon absorption, which makes the 6 g = 8 p peak in 
the susceptibility about an order of magnitude larger than the 6 g = 0 (Fig. 13) peak. 
For increasing values of 0 g , the two photon peak saturates very rapidly, since the 
single photon peak is approximately the same height as the two photon absorption. 
The one photon term also has a small gain peak as well as an absorption peak at the 
two photon gain peak (Fig. 14). 

Another difference between the cesium system and the previous two is the 
narrowness of the two photon gain peak. This is because in the cesium set of rates 
the polarization decay rates are due to spontaneous emission and not collisions. Since 
all of the rates in the plots are normalized to T x and the width is roughly 
determined by Tj/Tj, which is the relative rate of two-photon coherence decay. In 
the first two systems this number is unity, while in the cesium system it is 
Tj/Tj = 0.063. Fig. 15(a) shows how the 6 g = 6 p peak narrows as the two photon 
contribution becomes more dominant with increasing pump field detuning, while Figs. 
15(b), (c) shows how the width of the the peaks in the other two rate cases remains 
about the same. 

Figs. 15(a-c) also show that the absolute maximum of the susceptibility 
initially increases with pump detuning at constant pump field in both the cesium and 
the coupled rate FIR system, while there is a monotonic decrease in the equal rate 
case. Figs. 16-18 show the value of the absolute maximum of the susceptibilities as a 
function of increasing pump detuning for pump field Rabi frequencies of two, five, 
and ten in units of Tj 1 . In the equal rate case there is always a rapid monotonic 
decrease, and although the coupled rate FIR system rises to a peak, it then decreases 
rapidly with detuning. The cesium system has its maximum gain at a greater 
detuning than the FIR system and falls off much slower with detuning. This 
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suggests that the cesium may be used more effectively as a tunable laser medium 
than would be expected from looking at the equal rate case, since there is an 
appreciable gain peak over a significantly larger range of detunings. 

Fig. 19 shows the position of the peak signal frequency detuning as a 
function of pump detuning and the deviation from 8 g = 8 p as 8 g - 8 p is plotted vs. 
0 p . The peak location is the same for the three sets of rates, illustrating that for 

this range of pump power, the peak location depends only on 0 p and 8 p . At 8 p = 0 

the peak is located at 8. = 0„, the Rabi frequency. Only at very low powers where 
the Rabi splitting is not simply equal to the Rabi frequency do the rates affect the 
position of the peaks (Fig. 20). Fig. 21 shows the peak position vs. pump field 
Rabi frequency 0 p at 8 p = 0 and 8 p « 10. For 8 p = 0, the peak is at 8 g = 0 for 
low power, then splits at different powers for the three different rate cases, and 

eventually approaches 8 g = 0 p at higher powers. For 8 p = 10, the peak is at 8 g = 8 p 

and then as pump power is increased the peak shifts towards 8 g = 0 p . 
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III. Inhomogeneous Model 

In this section we will examine the effects of Doppler broadening in the 
cesium system for different values of the pump field detunings and Rabi frequencies. 
We will first see how the gain curves evolve as the detuning approaches and exceeds 
the Doppler width. Then, in the second half of this section, we will compare the 
predicted tuning curves for the Doppler broadened system vs. the homogeneous system 
for pump detunings and pump Rabi frequencies much greater than the Doppler 
width, as would be present in the doubled alexandrite pumped cesium vapor Raman 
laser. 

The Doppler broadened gain is given by the homogeneous response of a 
specific velocity group integrated over the Boltzman velocity distribution, which can 
be written as 


*Dop 



2 



-41n2(u/Au 

X h (8, T € u, 8 p - u)e 



du 


(26) 


where e * k g /k p , k p and k g are, respectively, the pump and signal field wave vectors, 
u is the angular frequency shift of the pump, and Au Dop is the full width at half 
max of the pump transition given by 
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The - and + are, respectively, the co- and counter-propagating beams. 

For cesium, the transition energies in wave numbers are 25709 cm' 1 for the 
bSj/jSPj/j transition and 1572 cm' 1 for the transition. Using these values 

and a temperature of 1000K.(^1 (the cesium vaporization temperature is 932K), Au Dop 
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and the wave vector ratio were calculated to be 

*>Dop = 54 - 5 T i X € = k s/ k p = °' 061 ( 28 > 

The integral was evaluated numerically. 

For small pump field and detunings^ within the Doppler width, the line 
shape exhibits a co- and counter-propagating asymmetry with the co-propagating peak 
located at c8 p and the counter-propagating peak located at -€8 p . The co-propagating 
peak is slightly higher and narrower than the counter-propagating peak by 
approximately 1 + e (Fig. 22), and as the pump Rabi frequency is increased above 
0 p - 7/5 the co-propagating peak splits, although the magnitude of the splitting 
is less than in the homogeneous case (Fig. 23). 

As the detuning is increased, the co-propagating split disappears, and a wide 
peak begins to appear at the Rabi shifted two photon resonance. The origin of this 
can be seen by approximating the homogeneous line as a sum of delta function 
responses at the one and two photon resonances. 

X h = V l6 (8 8 - A.) + V 2 8(8 s - 8 p - A 2 ) (29) 

where V x and V 2 are the areas of the two peaks, and A x and A 2 are the Rabi light 
shifts. Although the homogeneous line will only look like delta functions when 
(k f /k p ) Aw Dop »T., kj, this approximation shows some important features of the fully 
integrated line shape. 


-41n2 


Xd °p“ eAw, 


8 - A, 

S 1 


€ Aw, 


Dop 


-41n2 


Dop 


( 1 ¥e)Au Dop 


6 - 5 n‘ A 1 

s p 2 

Aw Dop 


^2 


( 30 ) 



-16- 


where the weak dependence of A x and A 2 on 8 p was neglected. The first term is just 
a gaussian near line center with the width of the laser transition eAu Dop . The second 
term near 8 g ~ 6 p is what the Raman peak becomes in the presence of Doppler 
broadening. It is a gaussian with the width of the pump transition ? the width of 
the laser transition, depending on co- or counter-propagating beams. Also, the peak 
value of the homogeneous two photon transition is diminished by l/(lT€)Au Dop , while 
the one photon peak is only diminished by l/€Au Dop . Thus, the 8 g = A x peak will be 
higher and narrower than the S g = 8 p + A 2 peak. Fig. 24 shows the fully integrated 
co- and counter-propagating gain curve for 8 p = lOOrj 1 and 0 p ** Tj 1 . It is 
interesting to contrast the Doppler integrated gain where the one photon peak is 
significantly higher than the two photon peak, with the homogeneous gain curve, in 
which the single photon peak is almost negligible in comparison to the two photon 
peak (Fig. 25). 

For higher pump powers the on resonance counter-propagating gain peak 
begins to split (Fig. 26), and as the pump Rabi frequency, B p , approaches the half 
width of the Doppler distribution, both the co- and counter-propagating gain curves 
split into two distinct peaks located at ±0 p . The evolution of the curves with 
increasing detuning proceeds much the same as the low pump power except that the 
co- and counter-propagating asymmetry persists on the single photon peak. 

We will now compare the predicted gain tuning curves of a doubled 
alexandrite pumped cesium vapor Raman laser to examine the effects of Doppler 
broadening. The alexandrite laser is a solid state laser, tunable from 7000A to 8000A. 
We will assume a pulsed laser with an energy of 150 mJ/pulse and a pulse width of 
T - 300 ns, focused to a spot size of 100 nm radius, which gives a peak Rabi 
frequency of 0 p - 3.6 x 1 0 5 rj 1 . The tuning range which will be examined is from 

25709 cm" 1 to 27200 cm" 1 , keeping in mind that the higher wave numbers will begin 
to pump the 9P^ 2 level, and thus will be most accurate closer to resonance. When 
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tuned to resonance at peak power, the susceptibility has two peaks Rabi split to 

approximately ±275 cm' 1 wave number detuning from the laser transition line center. 
The peak position is the same for both the homogeneous and inhomogeneous 

susceptibilities, although the peak gain value is much lower in the Doppler case than 
it would be if the system were not inhomogeneously broadened. As the pump is 
detuned the peak closest to the pump frequency begins to increase and move away 
from the center, while the other peak decreases and moves towards line center. Fig. 
27 shows the resonantly pumped susceptibility and the value of the peak as the 

detuning is increased. In Fig. 28 the position of the peak is given for varying pump 
detunings, with the diagonal line denoting the 8 g = 8 p line. The dotted line indicates 
the peak closer to line center, while the solid line indicates the "two photon" peak, 
which approaches 8 g = 8 p as the detuning becomes large compared with the Rabi 
frequency. 

When inhomogeneous broadening is taken into account, there is a drastic 
deviation from the homogeneous results. Although the peak position remains the 

same as in the homogeneous line, the peak which approaches 8 g = 8 p , widens and 
decreases in height from its 8 p - 0 value while the other peak narrows and increases 
sharply, even though the detuning and Rabi frequency are both huge in comparison 
to the Doppler width. Figure 29 shows how the peak changes with detuning from 
the S p - 0 value. 

This illustrates the necessity of including Doppler broadening in the gain 
calculation for such a system. The essence of the lower gain on the two photon 
transition is because of the near proportionality of peak position, 8 g , to pump 
detuning, 8 p . In contrast to the homogeneous calculations, including the calculation 
of [9,10,11] j n w hi c h each atom radiates at the same frequency, the Doppler 
broadened calculation includes a frequency shifted peak for each velocity group. This 
results in a two photon gain which is spread out with the width approximately equal 
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to the Doppler width of the pump transition and height proportionately diminished, 
as shown in Eq. (30). In addition, if the vapor cell is used as a one or two pass 
Raman amplifier, the line width of the output radiation will be that of the pump 
Doppler width, and not the homogeneous or power broadened linewidth. 
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IV. Conclusion 

We have examined the imaginary susceptibility on the laser transition 
(proportional to the gain or attenuation) of a three-level system, with arbitrary 
relaxation rates, detunings, and field strengths interacting with two laser fields, 
obtained analytically in the density matrix formalism. The solution assumed 
negligible spatial and temporal variation and monochromatic laser fields (laser line 
width much smaller than the homogeneous width). It was found that differing 

relaxation rates cause important qualitative differences, and some case, gain in one 
system may become absorption in another. 

The results were then applied to a cesium vapor laser pumped by a tunable 
pulsed doubled alexandrite laser. Since the pulse width is just slightly larger than 
the homogeneous line width, the system will be just on the limit of the adiabatic 
following regime. Therefore, the steady state results can only be loosely applied and 
must serve as a rough model for the position and relative strengths of the gain 

peaks. Another complication due to the pulsed pump arose because of the sensitivity 
of the gain peak position on pump field strength. This will result in "chirped" 

output pulses. However, since the peak gain for fixed detuning and pulsed pump 

field strength occurs at the peak of the pump pulse, the dominant contribution to the 
output gain will come from the peak field strength, resulting in a relatively narrow 
output pulse. 

The system was then Doppler integrated, since the Doppler width is many 
times the homogeneous width. It was found that even though the detuning and Rabi 
frequency were much greater than the Doppler width, it caused the gain to shift 
from the predominantly two photon (8„ - 8J homogeneous susceptibility result to a 
predominantly single photon contribution (8 S - 0). 
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Figure Captions 


Figure 1: Schematic of 3-level system. 


Figure 2: Energy levels of atomic cesium. 


Figure 3: x|j V T 1 M 2 € 0 is plotted versus 8 g T 1 with decreasing amplitude corresponding 

to increasing stimulated fields of 0 g T x = 0, 1, 2, 5 and (a) 0 T x = 1, 
Vl = 5, (b) 0 Tj = 5, SpTj = 5, (c) 0 p T x = 1, 8 p T 1 = 10, (d) 0 p T 1 = 5, 

Vi - 10 - 


Figure 4: (a) One and two photon terms with two photon term absorption near 

8 g T i " 0 of xJ! ft/T x M 2 € 0 and (b) x|j h/T 1 y 2 € 0 are plotted versus 8^ for 

8 p t x = 10 « Vi a 2 ’ 6 . T i = °- 


Figure 5: Same as Fig. 4 except for 8 t x = 10, 0 p T x = 2, 0 g T x = 1. 


Figure 6: Same as Fig. 4 except for 8 p T 1 = 10, 0 p T x = 2, 0 g T x = 5. 


Figure 7: Same as Fig. 4 except for 8 p T x = 10, 0 p T x = 2, 0 g T x = 10. 

Figure 8: (a) One photon term has higher amplitude of xJJ h/T 1 n 2 e Q and (b) 
Xjj fi/T 1 n 2 e 0 are plotted versus 8 g T x for 8 p T x = 0, 0 p r x = 0.2, 0 g r x = 0. 

Figure 9: (a) Single photon term is absorptive at 8 g = 0 and (b) xJJ fi/ T i n 2 e o are 
plotted versus 8 g T x for SpTj = 0, 0 p T 1 = 5, 0 g T x = 0. 

Figure 10: Same as Fig. 9 except for 8 p T 1 = 0, 0 p T x = 5, 0 g T x = 8. 

Figure 11: Same as Fig. 4 except with coupled FIR rates for 8 T, = 10, 0 T. = 2, 
*Ji = 0 - 


Figure 12: Same as Fig. 11 except for 8 p r 1 = 10, 0 p T x = 2, 0^ = 2. 


Figure 13: Same as Fig. 11 except with cesium rates for 8 T x = 

0 g T x = 0. 
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Figure 14: Same as Fig. 13 except for = 10, 0^ = 2, = 1. 

Figure 15: xJJ h/T 1 n 2 e 0 is plotted versus b a T l at values of 8^ = 2, 7, 10, 15 for (a) 

cesium rates and (b) equal coupled rates. 

Figure 16: Peak value of x£ fi/r 1 n 2 e o is plotted versus 6 at 0^ = 2, 0^ = 0 for 

(a) cesium rates, (b) equal coupled FIR rates, and (c) equal uncoupled 

rates. 

Figure 17: Peak value of xj h/T 1 n 2 e 0 is plotted versus 6^ at 0 p r x = 5, 0^ = 0 for 

(a) cesium rates, (b) equal coupled FIR rates, and (c) equal uncoupled 

rates. 

Figure 18: Peak value of x£ h/T 1 n 2 € 0 is plotted versus ITTj at 0 p T t = 10, 0^ = 0 

for (a) cesium rates, (b) equal coupled FIR rates, and (c) equal 
uncoupled rates. 


Figure 19: (a) peak position (8^) peak is plotted versus 8^ at 0 p T x = 10, 0 g T 1 = 0 

and (b) deviation from 8 g = 8 p is plotted versus 6 p T r 

Figure 20: (a) peak position (8^) peak is plotted versus 0 T 1 at 8^ = 0, 0 S T X = 0 

for (a) cesium rates, (b) equal coupled FIR rates, and (c), equal 
uncoupled rates. 

Figure 21: (a) peak position (8^) peak is plotted versus 0^ at S p T 1 = 10, 
0»Ti = 0. 


Figure 22: (a) Co-propagating and (b) counter-propagating values of x!J h/T 1 ti 2 e Q 

versus 8 g r 1 Doppler integrated at BpTj = 0.5, 0 g T x = 0, 8 p T 1 = 10. 

Figure 23: (a) Co-propagating and (b) counter-propagating values of fi/T l ^ 2 € o 

versus 8^ Doppler integrated at 0 Tj = 1.0, 0^ = 0, - 10. 

Figure 24: (a) Co-propagating and (b) counter-propagating values of x g fi/r,/x 2 € 0 

versus 6^ Doppler integrated at 0 p = 1.0, 0 S T 1 = 0, 8^ = 200. 


Figure 25: Homogeneous value of xJJ versus 8 g T 1 for 0 = 

8 p T i = 200. 
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Figure 26: (a) Co-propagating and (b) counter-propagating xJJ ft /f 1 M 2 e o versus 6 g T 1 

Doppler integrated at BpTj = 20, 0 g T x = 0, 8 p T x = 0. 

Figure 27: Vertical lines are a plot of the on resonant h/T l ii 2 € Q at 

0 g = 3.6 x lO 5 ^, and the curves are a plot of the maximum value of 
the respective homogeneous peaks versus 8 g in wave numbers as pump 
detuning is decreased below zero. 

Figure 28: A plot of the values of the two peaks of xJJ h/T l ii i e 0 at 0 g = 3.6 x 10 5 r 1 

versus 8„ in wave numbers. 


Figure 29: Vertical lines are a plot of the on resonant Doppler integrated 
Xjj ft/ t 1 m 2 « 0 at 0 g = 3.6 x 10 5 r r and the curves are a plot of the 
maximum value of the respective Doppler integrated peaks versus 8 g in 
wave numbers as pump detuning is decreased below zero. 


Table 1 


Transition Dipole Matrix Element (C-m) Spontaneous Lifetime (Sec] 


6s 1/2 6 Pl/2 

2.8x1 O' 29 

1.09 x 10-8 

6Si/ 2 7p 1 /2 

3.1 x 10- 30 

1.17 x 10' 7 

6s 1/2 8p 1/2 

1.4 x 10' 30 

3.59 x 10" 7 

6 Pi/2 7s 1/2 

2.6 x 10" 29 

4.54 x 10-8 

6 Pi/2 8s 1/2 

5.7 x 10-30 

1.68 x 10’ 7 

7s 1/2 7 Pl/2 

6.1 x 10- 29 

9.41 x 10-8 

7s 1/ 2 8pi/2 

6.6 x 10- 29 

7.30 x 10-9 

6Pl/2 8s i/2 

5.8x10- 29 

2.67x1 O' 7 

8 Si/ 2 8pi/2 

1.0x10-28 

2.77 x IQ' 7 

Effective Rates (sec -1 ) 

k, = 2.62 x 10 7 

V = 1.53 x 10 8 


k 2 — 3.06 x 10 7 

x^ 1 = 1 .43 x 10 8 


k 3 = 1.24 x 10 7 
k 4 as 2.60 x 1 0 7 

X3 1 =9.71 x 10 6 
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18536 



6s 1/2 


27632 9p 3/2 
27637 9p 1/2 


25792 8p 3/2 
25709 8p 1/2 


21947 7p 3/2 
21766 7p 1/2 


11732 6p 3/2 
11178 6p 1/2 
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